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ABSTRACT 

We discovered fine-scale structure within the scattering disk of PSR B0329+54 in observations with 
the RadioAstron ground-space radio interferometer. Here, we describe this phenomenon, characterize 
it with averages and correlation functions, and interpret it as the result of decorrelation of the impulse- 
response function of interstellar scattering between the widely-separated antennas. This instrument 
included the 10-m Space Radio Telescope, the 110-m Green Bank Telescope, the 14 x 25-m Westerbork 
Synthesis Radio Telescope, and the 64-m Kalyazin Radio Telescope. The observations were performed 
at 324 MHz, on baselines of up to 235,000 km in November 2012 and January 2014. In the delay 
domain, on long baselines the interferometric visibility consists of many discrete spikes within a limited 
range of delays. On short baselines it consists of a sharp spike surrounded by lower spikes. The average 
envelope of correlations of the visibility function show two exponential scales, with characteristic delays 
of Ti = 4.1 ± 0.3 /is and T 2 = 23 ± 3 /is, indicating the presence of two scales of scattering in the 
interstellar medium. These two scales are present in the pulse-broadening function. The longer scale 
contains 0.38 times the scattered power of the shorter one. We suggest that the longer tail arises from 
highly-scattered paths, possibly from anisotropic scattering or from substructure at large angles. 

Keywords: scattering - pulsars: individual B0329+54 — radio continuum: ISM — techniques: high 
angular resolution 


1. INTRODUCTION 

All radio signals from cosmic sources are distorted by 
the plasma turbulence in the interstellar medium (ISM). 
Understanding of this turbulence is therefore essential 
for the proper interpretation of astronomical radio ob¬ 
servations. The properties and characteristics of this 
turbulence can best be studied by observing point-like 
radio sources, where the results are not influenced by 
the extended structure of the source, but instead are di¬ 
rectly attributable to the effect of the ISM itself. Pulsars 
are such sources. Dispersion and scattering affect ra¬ 
dio emission from pulsars. Whereas dispersion in the 
plasma column introduces delays in arrival time that 
depend upon frequency and results in smearing of the 
pulse, scattering by density inhomogeneities causes an¬ 
gular broadening, pulse broadening, intensity modula¬ 
tion or scintillation, and distortion of radio spectra in 
the form of diffraction patterns. The scattering effects 
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have already been studied extensively theoretically (see 
e.g., Prokhorov et al. 19751 IRickett |1977 |Goodman fe 


Narayan|T9 89; Naraya n fe Goodman 1989| [STiishov et a T 

20031 and obs ervationally with ground VL B1 of SgrA* 

(ICwmn et al.||2014l) and pulsars (see, e.g., IBar tel et al. 

][ir - ' 


985| Desai et al. 1992 


Kondratiev et al. 


2007p, as well 


as with ground-space VLB1 of PSR B0329+54 ( Halca . 


Yangalov et al. 

2001 

) and the quasar 3C 273 ( RadioAs- 

iron, Johnson et al. 

2016 

). Whereas the VSOP pulsar 


1.7 GHz and with baselines of «25,000 km and less, 
ground-space VLBI with RadioAstron allows observa¬ 
tions at one-fifth the frequency, where propagation ef¬ 
fects are expected to be much stronger, and wit h base¬ 
lines ~10 times longer (Kardashev et al. 2013). Such 
observations can resolve the scatter-broadened image of 
a pulsar and reveal new inform ation about the scattering 
medium ( Smirnova et al.| [2014). 

In this paper, we study the scattered image of the pul¬ 
sar B0329+54 with RadioAstron. We demonstrate that 
the pulsar is detected on baselines that fully resolve the 
scattering disk. The interferometric visibility on these 
long baselines takes the form of random phase and am¬ 
plitude variations that vary randomly with observing fre¬ 
quency and time. In the Fourier-conjugate domain of 
delay and fringe rate, the visibility forms a localized, 
extended region around the origin, composed of many 
random spikes. We characterize the shape of this region 
using averages and correlation functions. We argue the¬ 
oretically that its extent in delay is given by the average 
envelope of the impulse-response function of interstellar 
scattering, sometimes called the pulse-broadening func¬ 
tion. We find that the observed distribution is well-fit by 
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Figure 1. Relations among the interferometric visibility V in var¬ 
ious domains, and functions derived from it. The fundamental ob¬ 
servable is the visibility in the domain of frequency v and time t , 
V(u,t)] this is known as the cross-power spectrum, or cross spec¬ 
trum. An inverse Fourier transform of v to delay r leads to the 
visibility V (r, t) ; this is the cross-correl ation function of electric 
fields in the time domain (see Equation |21| ). A forward Fourier 
transform of t to fringe rate f leads to V (r, /). A forward trans¬ 
form of t back to v produces V{y, /), and an inverse Fourier trans¬ 
form of / to t returns to V(u, t) The square modulus of V (r, t) is 
C(t , t). The cross-correlation function in r of Cr for right- and Cl 
for left-circular polarization is Krl{ At, t). We denote the Fourier 
transform by $, and quantities in the domain of frequency v by the 
accent ~. 


a model that is derived from an impulse-response func¬ 
tion that has two different exponential scales. We discuss 
possible origins of the two scales. 


2. THEORETICAL BACKGROUND 

Our fundamental observable is the interferometric visi¬ 
bility V. In the domain of frequency v, this is the product 
of electric fields at two antennas A and B : 


V AB (y,t) = E A (v,t)E%(v,t). 


(1) 


This representation of the visibility is known as the cross 
spectrum, or cross-power spectrum. Because electric 
fields at the antennas are complex and different, V A b 
is complex. Usually visibility is averaged over multi¬ 
ple accumulations of the spectrum, to reduce noise from 
background and the noiselike electric field of the source. 
The second argument t allows for the possibility that the 
visibility changes in time, as it does for a scintillating 
source, over times longer than the time to accumulate a 
single spectrum. Such a spectrum tha t changes i n time 
is known as a “dynamic spectrum” (Bracewell 2000). 
The correlator used to analyze our dat a, as discusse d in 
Sections [ 3 ] and [ 4 } calculates V A b(v, t) (Andrianov et al. 
2014|. Hureafter, we omit the baseline subscript Indicat- 
ing baseline AB in this paper, except in sections of the 
Appendix where the baseline is important. 

Under the assumptions that the source is pointlike, 
and that we can ignore background and source noise, 
the impulse-response function of interstellar scattering g 
determines the visibility of the source. A single delta- 
function impulse of electric Held at the source is received 


as a function g(t e ) of time t e at the observer. Here, t e 
is Fourier-conjugate to v and varies at the Nyquist rate. 
The visibility is the product of Fourier transforms of g 
at the two antennas: 

V AB = g A g* B ( 2 ) 


where g is the Fourier transform of g{t e ). 

We denote the typical duration of g(t e ) as r sc , the 
broadening time for a sharp pulse. Within this time span, 
g{t e ) has a complicated amplitude and phase. The func¬ 
tion g(t e ) changes over longer times, as the line of sight 
shifts with motions of source, observer, and medium. 
This change takes place on a timescale t sc , and over a 
spatial scale S sc ■ The shorter and longer timescales t sc 
and t sc lead to our use of dual time variables: t e , of up 
to a few times r sc and Fourier-conjugate to and t, of a 
fraction of t sc or more and Fourier-conjugate to /. This 
duality is commonl y ex pressed via the “dynamic spec¬ 
trum” (see Section A.2|. If the scattering material re¬ 
mains nearly at rest while the line of sight travels through 
it at velocity Vj_, then one spatial dimension in the ob¬ 
server plane maps into time, and 


tsc — ( S S c / 1 L 


( 3 ) 


The averaged square modulus of g is the pulse¬ 
broadening function G = (g(t e )g(t e )*)s- Here, the sub¬ 
scripted angular brackets (...)§ indicate an average over 
realizations of the scattering. This function is the aver¬ 
age observed intensity for a single sharp pulse emitted 
at the source. An average over time is usually assumed 
to approximate the desired average over an ensemble of 
statistically-identical realizations of scattering. 

We derive a number of representations of the visibil¬ 
ity and quantities derived from it, and show that these 
provide straightforward means to extract the impulse- 
response function. These functions are summarized in 
Figure lj and discussed briefly here, and in detail in Sec¬ 
tion [A of the Appendix. In particular, visibility in the 
domain of delay r and time t is V{r,t). This is the cor¬ 
relation function of e lectric field at the two antennas A 
and B (Equation [2l]), and is the inverse Fourier trans¬ 


form of V{v, t) from v to r. We are also c oncern ed with 
the square modulus of V(r,t) (see Section A.3.2): 


C(r,t) = |U(r,t)| 2 (4) 


We calculate C for right- and left-circular polarizations 
separately, and then correlate them in delay r to form 
Krl, the cross-correlation between polarizations: 

Kr L (At, t) = ^Y jt C r {tA)C l {t + Ar,t) (5) 

Here, Krl is the correlation of a single measurement of 
Cr and Cl, and N is the number of samples in Cr and 

C L . 

When averaged over many realizations of the scat¬ 
tering material, ( Krl)s is related to the statistics 
of the pulse-broadening function G. Most commonly, 
the average over many realizations of scattering ma¬ 
terial is approximated by averaging over a time much 
longer than f sc ; for this reason we omit the time ar¬ 
gument for (. Krl(At))s ■ Equivalently, evaluation of 
{Krl{At , frnax)) ~ (Cr[t, fmax)C l{t T At, i /* maa ,)) T at 
the fringe rate f m ax of the maximum magnitude of Krl 
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Table 1 

Diary of observations 


Epoch of 

Time 

Ground 

Polarizations 

Scan 

Observations 

Span 

Telescopes 


Length 

2012 Nov 26 through 29 

1 hr/day 

GB 

RCP+LCP 

570 s 

2014 Jan 1 and 2 

12 hr 

WB, KL 

RCP 

1170 s 


Table 2 

Observations on Earth-Space Baselines 


Epoch 

Projected 
Baseline Length 
(10 3 km) 

RA 

Observing Time 
(minutes) 

2012 Nov 26 

60 

60 

2012 Nov 27 

90 

60 

2012 Nov 28 

175 

60 

2012 Nov 29 

235 

60 

2014 Jan 1 

20 

60 

2014 Jan 2 

70 

100 

2014 Jan 2 

90 

120 


yields the same time average. For this theoretical discus¬ 
sion, fmax = 0; for practical observations, instrumental 
factors can offset the fringe rate from zero, so that f max 
provides the most reliable time average. 

For a baseline that extends m uch further than the scale 
of scattering S sc (see Equation [32|): 


{Krl{t))s = G(i~) 0 G-(t) 
+ (l if t = 0 ) 


G(t)®G_(t) (6) 


Here, we introduce the symbol <S> to indicate convolution, 
and denote the time-reverse of G as G_(r) = G(— t). 

Our analy sis method differs somewhat from Smirnova 
et al. (2014), who used structure functions of intensity, 
visibility, and visibility squared to study scattering of 
pulsar B0950+08 on an extremely long baseline to Ra¬ 
dioAstron. The two methods are closely related theoreti¬ 
cally. Structure functions are particularly valuable when 
the characteristic bandwidth approaches the instrumen¬ 
tal bandwidth, and can be extended to cases where the 
signal-to-noise ratio is low, as they discuss. 

3. OBSERVATIONS 

The observations were made in two sessions: the first 
for one hour each on the four successive days November 
26 to 29, 2012, and the second for a total of 12 hours 
on the two days January 1 and 2, 2014. The first ses¬ 
sion used the 10-m RadioAstron Space Radio Telescope 
(RA) together with the 110-m Robert C. Byrd Green 
Bank Telescope (GB). The second session used the RA 
together with the 14 x 25-m Westerbork Synthesis Ra¬ 
dio Telescope (WB), and the 64-m Kalyazin Radio Tele¬ 
scopes (KL). Both right (RCP) and left circular polariza¬ 
tions (LCP) were recorded in November 2012, and only 
one polarization channel (RCP) was recorded in January 
2014. Because of an RA peculiarity at 324 MHz, the 
316-332 MHz observing band was recorded as a single 
upper sideband, with one-bit digitization at the RA and 
with two-bit digitization at the GB, WB, and KL. Science 
data from the RA were transmitted in real time to the 


telemetry station in Pushchino (Kardashev et al. 2013) 
and then recorded with the RadioAstron data recorder 


(RDR). This type of recorder was also used at the KL, 
while the Mk5B recording system was used at the GB 
and WB. Table [T] summarizes the observations. 

The data were transferred via internet to the Astro 
Space Center (ASC) in Moscow and then processed with 
the ASC correlator with gating and dedispersion applied 
(Andrianov et al. 2014). To determine the phase of the 
gate in the pulsar period, the average pulse profile was 
computed for every station by integrating the autocor¬ 
relation spectra obtained from the ASC correlator. The 
autocorrelation spectra Vaa(+ t) are the square modulus 
of electric field at a single antenna. 

In November 2012 the projected baselines to the space 
radio telescope were about 60, 90, 175, and 235 thou¬ 
sand kilometers for the four consecutive days, respec¬ 
tively. Data were recorded in 570-second scans, with 
30-second gaps between scans. In January 2014 the pro¬ 
jected baselines were about 20, 70, and 90 thousand kilo¬ 
meters during the 12-hour session. Data were recorded in 
1170-second scans. The RA operated only during three 
sets of scans of 60, 100 and 120 min each, with large 
gaps in between caused by thermal constraints on the 
spacecraft. The auto-level (AGC), phase cal, and noise 
diode were turned off during our observations to avoid 
interference with pulses from the pulsar. Table [2] gives 
parameters of the Earth-space baselines observed - ^ 


4. DATA REDUCTION 
4.1. Correlation 

All of the recorded data were correlated with the ASC 
correlator using 4096 channels for the November 2012 
session and 2048 channels for the January 2014 session, 
with gating and dedispersion activated. The ON-pulse 
window was centered on the main component of the av¬ 
erage profile, with a width of 5 ms in the November 2012 
session and 8 ms in the January 2014 session. These 
compare wit h a 7 -ms pulse width at 50% of the peak 
flux density (Lorimer et al. 1995). The OFF-pulse win¬ 
dow was offset from the mam pulse by half a period and 
had the same width as the ON-pulse window. The cor¬ 
relator output was always sampled synchronously with 
the pulsar period of 0.714 s (single pulse mode). We 
used ephemerides co mputed with the prog ram TEMPO 
for the Earth center (Edwards et al. 2006). The results 
of the correlation were tabulated as cross power spectra, 
V{y,t) 1 written in standard FITS format. 


4.2. Single-Dish Data Reduction 

Using autocorrelation spectra at GB, KL, and WB, 
we measured the scintillation time t sc and bandwidth 
A v sc = 1/27 tt sc . The results are given in Table [3] Our 
analysis using interferometric data, for which the noise 
baseline is absent and the spectral resolution was higher, 
is more accurate for the constants T\ and 12 as discussed 
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Table 3 

Measured Scattering Parameters of PSR B0329+54 


Epoch 

^sc 

Ak'sc 

®nr 

w nf 

n = 1/fci 

r 2 = 1 /k 2 


(s) 

(kHz) 

(ns) 

(mHz) 

(rs) 

(rs) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

Nov 2012 

114 ±2 

15 ±2 

50 ± 5 

20 ± 2 

4.1 ± 0.3 

23 ±3 

Jan 2014 

102 ±2 

7± 2 

43 ±3 

25 ±3 

7.5 ± 0.3 

- 


Note. — Columns are as follows: (1) Date of observations, (2) 
Scintillation time from autocorrelation spectra as the half width at 
1/e of maximum, (3) Scintillation bandwidth from single-dish au¬ 
tocorrelation spectra as the half-width at half maximum (HWHM), 
(4) HWHM of a sine function fit to the central spike of the visibility 
distribution along the delay axis, (5) HWHM of a sine function fit 
to the central spike of the visibility distribution along the fringe 
rate axis, (6) Scale of the narrow component of |A'rl(At)| (Section 


5.2.31, (7) Scale of the broad component of |ATrr(At)| (Section 


below, so we quote those values in Table [3] 


4.3. VLBI Data Reduction 
The ASC correlator calculates the cross-power 


trum, V(u,t), as discussed in Sections [ 2 ] and A.3.1 


spec- 
The 

resolution of the resulting cross-power spectra is 3.906 
kHz for the 2012 observations and 7.812 kHz for the 2014 
observations. Because the scintillation bandwidth was 
comparable to the channel bandwidth for the 2014 ob¬ 
servations, as shown in Table [3j and because the single 
recorded polarization at that epoch prevented us from 
correlat ing po larizations to form A/rr, as discussed in 
Section 5.2.3[ we focus our analysis and interpretation 
on the 2012 observations. 


5. ANALYSIS OF INTERFEROMETRIC VISIBILITY 

We investigated the scattering of the pulsar from the 
visibility in the delay-fringe-rate domain, V(t, /). We 
studied the statistics of visibility V (r, /) as a function 
of delay, fringe rate, and baseline length. If there were 
no scattering material between the pulsar and the ob¬ 
server, we would expect for | V (r, /) | one spike at zero 
delay and fringe rate with magnitude that remains con¬ 
stant as a function of baseline length, and with width 
equal to the inverse of the observed bandwidth in delay, 
and the inverse of the scan length in fringe rate. Scat¬ 
tering material in between changes this picture. First we 
expect the spike at zero delay and fringe rate to decrease 
in magnitude with increasing baseline length, perhaps 
to the point where it would become invisible. Second, 
we expect additional spikes to appear around the spike 
at zero delay and fringe rate. The distribution of these 
spikes give us invaluable information about the statistics 
of the scattering material. 

As we discuss in this section, we fitted models to the 
distribution of visibility, as measured by the correlation 
function ATrr, and thus derived scintillation parameters 
that describe the impulse-response function for propaga¬ 
tion along the line of sight from the pulsar. We also com¬ 
puted the maximum visibility as a function of projected 
baseline length, as we discuss in detail in a separate pa¬ 
per (Paper II: Popov et al., in preparation). 

For strong single pulses the visibility in the cross spec¬ 
tra, V(u,t), had signal-to-noise ratios sufficiently large 


for a useful analysis. However, we decided to ana¬ 
lyze the data from the time series of multiple pulses. 
Fourier transform of the cross spectrum, V(iy,t), to the 
delay/fringe-rate domain yields V(u, f ) and concentrates 
the signal into a central region, and thus provides a high 
signal-to-noise ratio. The sampling rate of individual 
cross spectra in the time se ries was the pulse period of 
0.714 s, as noted in Section |4~l) The time span of cross 
spectra used to form V(r, /) varied, ranging from 71.4 s 
to 570 s, depending upon the application. 

5.1. Distribution of visibility 

In Figure[2]we display the magnitude of the visibility in 
the delay/fringe-rate domain, |V(t, /)|, for a 500-s time 
span. The data were obtained on 29 November 2012 in 
the RCP channel for a projected 200MA GB-RA base¬ 
line. The cross spectra, V(i/, t), from which we obtained 
| V(t, /) | were sampled with 4096 spectral channels across 
the 16-MHz band, at the pulsar period of 0.714 s; con¬ 
sequently, the resolution was 0.03125 /is in delay, and 
2 mHz in fringe rate. As Figure [2] shows, no dominant 
central spike is visible at zero delay and fringe rate, as 
would be expected for an unresolved source. Our long 
baseline interferometer completely resolves the scatter¬ 
ing disk. Instead we see a distribution of spikes around 
zero delay and fringe rate that is concentrated in a rela¬ 
tively limited region of the delay-fringe rate domain. The 
locations of the various spikes appear to be random. Be¬ 
cause the scattering disk is completely resolved on our 
long baseline, we conclude that the spikes are a conse¬ 
quence of random reinforcement or cancellation of paths 
to the different locations of the two telescopes, and hence 
interferometer phase. 

In Figure[2] the distribution of the magnitude of visibil¬ 
ity is relatively broad along the delay axis and relatively 
narrow along the fringe rate axis. The extent is limited in 
delay to about the inverse of the scintillation bandwidth, 
t S c = 1/27tAr sc ; and in fringe rate to about the inverse 
of the diffractive timescale tdis- Within this region, the 
visibility shows many narrow, discrete spikes. If statis¬ 
tics of the random phase and amplitude of scintillation 
are Gaussian, and the phases of the Fourier transform 
randomize the different sums that comprise the visibility 
in the delay-fringe rate domain, then the square modulus 
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Figure 2. Magnitude of visibility in the delay-fringe rate domain |V(r, /)|, for a 500-s time span on 29 November 2012 in the RCP 
channel, on the RA-GB baseline. Visibility is normalized for autocorrelation: |V(0,0)| = 1. The axes show instrumental offsets, including 
about 6 fis in delay. Top: three-dimentional representation; bottom: two-dimentional representation. 
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Figure 3. Examples of the fine structure of the magnitude of vis¬ 
ibility, | V(t, /max)I, as a function of delay t, with fringe rate fixed 
at the maximum of the delay-fringe rate visibility near zero fringe 
rate, fmax • From lowermost to uppermost, the curves correspond 
to progressively longer baselines, with the telescopes indicated and 
the approximate baseline projections given in MA in parentheses. 
Curves are offset vertically, and the upper 2 magnified as the ver¬ 
tical scale indicates, for ease of viewing. All curves show 71.4 s 
of data. Uppermost curve is from 2012 November 29; the rest 
are from 2014 January, when multiple ground telescopes provided 
shorter baselines. Note variation in scattering time between epochs 
as given in Table [3] The uppermost panel is the cross-section of 
the data shown invigu re [2l but for 71.4 s integration. Visibility is 
normalized as in Figure|2| t he best estimate of instrumental delay 
has been removed for each curve. 


of V (t, /) should be drawn from an exponential distribu¬ 
tion, multiplied by the envelope defined by the determin¬ 
istic part of the impulse-response function, as discussed 
in the Appendix. 

Along the delay axis, \V(r,f)\ takes the general form 
suggested by Figures [2] and [3] a narrow spike surrounded 
by a broad distribution. We found that the central spike 
takes the form of a sine function in both delay and fringe 
rate coordinates, as expected for uniform v isibility across 
a square passband (Thompson et al. 2007). The widths 
are somewhat larger than values expected from observ¬ 
ing bandwidth of 16 MHz and time span of 71.4 s, of 
w nT = 31.25 ns and w n f = 14 mHz respectively, proba¬ 
bly because of the non-uniformity of receiver bandpasses 
and pulse-to-pulse intensity variations, respectively. The 
broader part of the distribution takes an exponential 
form along the fringe-rate axis in this case; more gen¬ 
erally, the form can be complicated, particularly over 
times longer than 600 s. Traveling ionospheric distur¬ 
bances may affect the time behavior of our 92-cm obser¬ 
vations; in particular, they may be responsible for the 
20 to 25 mHz width of the narrow component in fringe 
rate, as noted in Table[3] We do not analyze the broader 
distribution in fringe rate further in this paper; we will 
discuss this distribution, and the influence of traveling 
ionospheric disturbances, in a separate publication (Pa¬ 
per III, Popov et al. in preparation). Because of the rel¬ 
atively small optical path length of the ionosphere, even 
at A = 9 2 cm , they cannot affect the cross spectrum 
( Hagfors||1976 ). 

The distribution of the magnitude of the visibility in 
delay/fringe-rate domain changes with baseline length. 
Figure [3] displays cross-sections through the maximum 


of the distribution of magnitude for a range of baseline 
lengths, as a function of delay. The maxima lie near 
zero fringe rate, as expected. Under the plausible and 
usual assumption that the correct fringe rate lies at the 
fringe rate, / max , where the distribution peaks, the cross- 
section represents the visibility averaged over the time 
span of the sample: 


V(T,f nm ) = (V{T,t)) t 


( 7 ) 


The top panel of Figure [3] shows this cross-section 
through Figure [2] The next lower panel shows the cross- 
section for the slightly shorter KL-RA baseline. The 
three lower plots give the equivalent cross-sections for 10 
times and 100 times shorter projected baselines. These 
three short-baseline cross-sections are qualitatively dif¬ 
ferent from the long baseline cross sections: the visibility 
has a central spike resulting from the component of the 
cross-spectrum that has a constant phase over frequency, 
as well as the broad distribution from the component 
that has a varying phase over frequency. The central 
spike is strongest for the shortest baseline and weaker 
for the next longe r basel ines, as ex pected based on the 
results of Sections |A.3.1| and |A.3.2| At very long base¬ 
lines the central spike is absent even after averaging the 
visibility over the whole observing period, and only the 
broad component is present. As expected from Figure 
[2j in the delay/time domain the broad component ap¬ 
pears as spikes distributed over a range of about 10 /us 
in delay. These spikes keep their position in delay for 
the scintillation time of about 100 to 115 s, as listed in 
Tabled 

The character of the broad component changes with 
baseline length as well: mean and mean square visibility 
are the same for short and long baselines; but excursions 
to small and lar ge visib ilities are more common for a long 
baseline (Gwinn 2001 Eq. 12). 


5.2. Averages and Correlation Functions 

Averages of the visibility, and averages of the correla¬ 
tion function of visibility, extract the parameters of the 
broad and narrow components of visibility. Such aver¬ 
ages approximate the statistical averages discussed in 
Sections [2] and [A] They seek to reduce noise from the 
observing system and emission of the source, as well as 
variations from the finite number of scintillations sam¬ 
pled, while preserving the statistics of scintillation. The 
averages and correlation functions allow the inference of 
parameters of the impulse-response function of propaga¬ 
tion from the statistics of visibility. 


5.2.1. Square Modulus of Visibility C 

The mean square modulus of visibility, ( C(t))s = 
(|U(t)| 2 )s, provides useful and simple characterization 
of visibility. To approximate the average over realiza¬ 
tions of scattering (...)§, we average over many samples 
in time t and over bins in delay r. We realize the av¬ 
erage over time by evaluating V(r,f) at the fringe rate 
of maximum amplitude f m ax , as discussed in Section [2] 
We also average over 16 lags in delay r. The result^ 
ing average shows a broad component surrounding the 
origin; on shorter baselines, it shows a spike at the ori¬ 
gin. The broad distribution samples the properties of 
the fine structure seen in Figures [2] and [3j and the spike 
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Figure 4. Cross-section of the mean square visibility in the 
delay/fringe-rate domain ( C(r))s = | V(r, fmax)\ 2 along the de¬ 
lay axis, at the fringe rate fmax where the magnitude of visibility 
peaks, close to zero mHz. The visibilities for the GB-RA baseline 
on 2012 Nov 28 at 21:40 UT are shown as open circles. The visibil¬ 
ities were computed by an inverse Fourier transform of the spectra, 
V(v,t), over 71.4 s time spans, and then by averaging over 6 ob¬ 
serving scans, each 570 s long. They were then further averaged in 
delay, over 16 points or 0.5 /xs, to smooth fluctuations. The dashed 
horizontal line shows the offset contributed by background noise. 
The sol id gray line shows the reconstructed form given by Equa¬ 
tion ^! offset by the noise level, with parameters taken from the fit 
shown in Figure [5] The light dashed curve shows only the narrow 
component of the two-exponential model. Units of visibility are 
correlator units. 

to those seen on the shorter baselines in Figure [3] We 
argue in Section |Al that the spike in (C(r))s is related 
to the average visibility, and the broad component to the 
impulse-response function. 

Figure [4] shows an example of the broad component of 
( C(t))s■ This is estimated as \V(t, fmax)\ 2 , by selecting 
the peak fringe rate fmax to average in time for each of 6 
scans, averaging the results for the scans, and averaging 
over 16 lags of delay to smooth the data. These averag¬ 
ing procedures serve to approximate the average over an 
ensemble of realizations of scattering. Background noise 
adds complex, zero-mean noise to V(t, /), with uniform 
variance at all lags; this adds a constant offset to the 
average (C(t)) s = |V(r, f ma x)\ 2 - 

5.2.2. Correlation Function K 

Using Equation[5j we estimated {Krl(At))s, the aver¬ 
aged cross-correlation function between the square mod¬ 
ulus of right-circular polarized (RCP) and of left-circular 
polarized (LCP) of visibility in the delay domain. (Note 
that (Krl(At))s is not the correlation function of the 
average ( C)s , but rather the average of the correlation 
function (C R <g) Cl)s ■) Because the background noise 
in the two circular polarizations is uncorrelated, they do 
not contribute an offset to {Krl(At))s- This allows us 
to follow the effects of the impulse-response function to 
much lower levels than for ( C)s■ The correlation func¬ 
tion (Krl(At))s is thus less subject to effects of noise, 
and is more sensitive to the broad component of the dis¬ 
tribution, than ( C) s■ 

To compute an estimate of ( Krl(At))s, we calculated 
the squared sum of real and imaginary components of 
U(r, f), the inverse Fourier transform of the cross-power 
spectrum. We formed these for each strong pulse, and 



Figure 5. An example of the correlation function (A'(Ar, / ma x))t 
on 2012 Nov 28, averaged over 570 s starting at 21:40:00 UT. The 
data were normalized by the square root of Krr and Krr at 
Ar = 0. The bes t-fit ting parameters for a 2-exponential fit of the 
form of Equation |35| are as indicated. 

normalized them by the autocorrelation functions at each 
antenna. From these we formed the un-averaged cor¬ 
relation function K R l(At, t). We then averaged Krl 
over 570-sec scans to form (A^i(Ar)) t . Averaging in 
the time domain approximated an average over realiza¬ 
tions of the impulse-response function for the scattering 
medium. Each 570-sec scan included 100 to 250 strong 
pulses, yielding one averaged sample of (K rl (At, t)) t for 
each scan. We obtained 22 measurements in total, with 
6 samples of (Arl(At, t)) t for November 26 , 28, 29 ob¬ 
serving sessions. We obtained only 4 such samples for 
November 27 because of no significant detections of V 
for two scans on that date. 

5.2.3. Two Exponential Scales 

Examination of the averaged cross correlation function, 
(K R L{AT,t)) t , revealed a spike at the origin and two 
exponential scales for the broad component, a large one 
and a small one. Figure [5] shows an example. 

The spike at the origin arises from the fine structure of 
scintillation in the broad component of visibility, as seen 
in Figures [2] and [3] This structure is identical in right- 
and left-circular polarizations, so its correlation leads to 
the spike. 

The two exponential scales are apparent as the slopes 
of the steeper and narrower parts of the distribution. 
We see these two scales even for single pulses, which 
are strong enough to show the two-scale structure. We 
did not observe these scales without doubt in spectra 
from single-dish observations, because the resulting cor¬ 
relation functions are more subject to noise, gain fluc¬ 
tuations, and interference. The scales are both present 
for ( C)s , but the longer scale is seen more clearly in 
(Krl{At))s (as comparison of Equations [34] through [36] 
shows). 

5.2.4. Model Fit 

We formalized the two exponential scales seen for 
(K(At))s with a model fit. The model assumed a pulse¬ 
broadening function with two exponential scales. Under 
this assumption, a short pulse appears at the observer 
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with average shape: 


(Aik\e klT + A 2 k 2 e fc2T , r> 0 
\0, r < 0 


( 8 ) 


The pulse rises rapidly, and falls as the sum of the two 
exponentials. 

The assumed form for G leads to predictions for the 
forms of (C) and ( K ), as discussed in Section [B] For 
(C), we expect a cusp at the origin, and two exponen¬ 
tials with scales ki and k 2 and different weights on either 
side. For ( K ), correlation smooths the cusp at the ori¬ 
gin, producing a smooth peak, with the same exponential 
scales appearing to either side. 

Figure 5] shows the best-fitting model of this form for 
the data shown there. This model has parameters: 


Ai/A 2 = 0.33 (9) 

ki = 1/4.3 p s (10) 

k 2 = 1/23 ps. (11) 

The model reproduces the two scales, and the smooth 
peak, well. The model also predicts the magnitude of the 
spike accurately, with zero average visibility pab = 0. 

The model shown in Figure ijshows the model for C, 
reconstructed using Equation 34] with parameters from 
the fit to Figure [5] The two scales appear in the model, 
although the offset from noise contributes at large delay 
r. As the figure shows, a single exponential does not fit 
the model well: the narrow component is satisfactory at 
small r, but falls well under the data at larger r. A high¬ 
winged function such as a Lorentzian can fit C well, but 
the rounded peak leads to a very wide peak for K that 
cannot match the data, and the inversion to a G(t) that 
remains finite, and is zero for t < 0 as causality demands, 
is problematic. 

The best-fitting scales and the magnitudes of the two 
contributions varied from scan to scan, but in a manner 
that was consistent with our finite sample of the scin¬ 
tillation pattern, and the inhomogeneous averaging of 
pulses with different intensities. We show a histogram 
of the results of our fits to 570-sec intervals in Figure 
[6] On 26 to 29 November 2012, the shorter scale av¬ 
eraged to ti = 4.1 ± 0.3 ps, and the longer scale to 
t 2 = 22.5 ± 2.9 ps. The scales had a relative power 
of A 2 /Ai = 0.38. 


6 . DISCUSSION 

On a long baseline that fully resolves the scattering 
disk, as Figure [3] shows, we observe multiple sharp spikes 
in the visibility V (t, /) is a consequence of the variation 
of the amplitude and phase of visibility. (See also Paper 
II, Popov et al. in preparation). The characteristic re¬ 
gion of that variation, At- A/, reflects the product of the 
inverses of the scintillation bandwidth At ss 1/2-k/S.v sc 
and the scintillation timescale A / « l/27rt sc . These 
quantities are the width in time of the impulse-response 
function, and the time for the impulse-response func¬ 
tion to change as the line of sight to the observer moves 
through the scattering material. 

Detailed examination of the correlation function of vis¬ 
ibility A'^jt(At, t) reveals the presence of two character¬ 
istic, exponential scales. Both scales are visible in the 
single-pulse correlation functions of right and left circu- 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 

Amplitude 

Figure 6. Upper panel: The distribution of long and short time 
scales for exponential scales of (lf^(Ar,/ max ))t. Each pair of 
scales was measured for a 570-s interval on one of the 4 consec¬ 
utive observing days in 2012. Lower panel: The distribution of 
magnitudes of long and short time scales. 


lar polarization, as well as in the correlation function 
(K RL (/S.T)) t averaged over 570 s shown in Figure [5] For 
an assumed screen distance of half the p ulsar distance of 


1.03iol2 kP c (|Brisken et al. 2002), the two scales 


cm 


cm 


( 12 ) 


D = 

correspond to diffractive scales of: 


tdi = £-\ — = 2.3 xlO 9 

Z7T V CTi 

(■d2 = 7T\f^- = !-0 x 10 9 

27T V ct 2 

The diffractive scale is the lateral distance at the screen 


where phases decorrelate by a radian (Narayan 1992). 
The refractive scale gives the scale of the scattering disk: 


cm 


l r i = i/ctiD = 1.9 x 10 
l r2 = \Jct 2 D = 4.6 x 10 13 cm 


(13) 


In contrast, Britton et al. (19981 measured angular 


broadening for PSR B0329+54' of 9h < 3.4 mas at v = 
325 MHz, where 9h is the full width of the scattered im¬ 
age at half the maxium intensity. This corresponds to a 
refractive scale of £ r = ( 0h/V& ln2)D/2 < 1.1 x 10 13 cm. 
This upper limit is somewhat smaller than the values ob¬ 
tained from our observations, even if one takes into ac¬ 
count the facts that the larger scale contains only 0.38 
of the power of the shorter one, and that the scattering 
material may be somewhat closer to the pulsar than to 
the observer. 


6.1. Previous Observations 


Shishov et al. (2003) studied the scattering properties 
of FS’R" B"03'2’9+54 in detail, using single-antenna obser¬ 
vations at 102 MHz, 610 MHz, 5 GHz, and 10.6 GHz 
to form structure functions of the scintillation in time 
and frequency on a wide range of scales. They con- 
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eluded that the scattering material has a power-law spa¬ 
tial spectrum with index a + 2 = 3.50 ± 0.05, marginally 
consistent with the value of 11/3 expected for a Kol¬ 
mogorov spectrum, with an ou ter scale o f 2 x 10 11 m < 


L 0 < 10 17 m. Using VLBI, Bartel et al. (1985) ob¬ 
served PSR B0329+54 at 2.3 tJHz and set limits on 
the separation on the emission regions corresp ondin g 
to dif feren t components of the pulse profile. |Yangalov 
et al. (20011 observed PSR B0329+54 at 1.650 GHz with 
ground-space baselines to HALCA, and found that the 
source varied strongly with time. They ascribed this vari¬ 
ation to scintillation, with the scintillation bandwidth 
comparable to the observing bandwidth at their observ¬ 
ing frequency. Self-calibration with timespans less than 
the scint illation time returned a pointlike image, as ex¬ 
pected. Semenkov et al. (2004) analyzed these data, 
including ground-ground baselines. They studied both 
single-antenna autocorrelation functions Vaa{At) and 
cross-correlation functions Var(At). They detected two 
timescales for the scintillation pattern, of 20 min and 1 
min. They found that the properties of scattering could 
not be explained by a single, thin screen, and further that 
velocities indicated relative motions within the scattering 
medium. 


Popov & Soglasnov (1984) had previously observed two 
coexisting scales ot scattering for PSR B0329+54 . They 
found scintillation bandwidths of Aiq = 115 Hz and 
Ab 2 = 750 Hz, measured as the 1/e point of the cor¬ 
relation function of intensity at an observing frequency 
of 102 MHz, using the Large Cophase Array of Puschino 
Observatory. The ratio of these scales, Ai^/Aiq = 6.5 
is larger than the ratio of fci/fc 2 = 5.5 that we observe. 
Scaled to our observing frequency of 324 MHz, using the 
Av oc b 22 / 5 scaling appropriate for a Kolmogorov spec¬ 
trum, and converting from Ai/ to r using the uncertainty 
relation r = 1 / 27 tAu, we find that these values corre¬ 
spond to 1.3 and 8 /zs, respectively, about a factor of 3 
smaller that the scales we observe. Of course, interpo¬ 
lation over a factor of 3 in observing frequency and the 
different observing techniques may introduce biases, and 
scattering parameters likely vary over the years between 
the two measurements. Two scales of scattering have 


also been observed fo r other pulsars (Gwinn et al. 
Smirnova et al.|[2014 ). 


2006 


6.2. Origin of Two Scales 

Two scales of scattering may be a consequence of a 
variety of factors. Non-Gaussian statistics of scattering 
can produce multiple scales, although this usually ap¬ 
pears as a continuum of scales rather than two different 
individual sc ales, as in a power-la w distribution or a Levy 
flight (Boldyrev & Gwinn 20031. A Kolmogorov model 
for scattering in a thin screen does not fit as well as our 
model based upon a two-exponential impulse-response 
function, or even as one based upon one exponential. A 
model with two discrete scales appears to fit our data 
better. 

One explanation is anisotropic scattering. This can 
produce two scales, corresponding to the major and mi¬ 
nor axes of the scattering disk, as discussed in Section [C] 
of the Appendix. The ratio of the scales of fc 2 /fci = 5.5 
corresponds to the parameter a 2 = 57, and an axial ra¬ 
tio of 02/di = Vl + a 2 = y / 2(fc 2 /fci) 2 — 3 = 7.4. In a 


simple model for anisotropic scattering in a thin screen, 
we expect the ratio of power in the scales to be approxi¬ 
mately \J 1 + a 2 /2-k{2 + a 2 ) ~ 0.40, as shown in Section 
O of the Appendix. This compares well with our ob¬ 
served ratio of A 1 /A 2 = 0.38. However, our observations 
for ground-space baselines at a variety of orientations do 
not show anisotropy. A variety of models, involving ma¬ 
terial with varying anisotropy distributed along the line 
of sight, and strong anisotropy that slips between our 
long baselines, might match our data. 

A second explanation is the complicated structure ob¬ 
served within dynamic spectra: most common ly observed 
as “scintillation arcs” (Stincbring et al. 2001). Recently, 
it has been suggested that this structure arises from in¬ 
terference among subimages, resulting fr om refract ion by 
interstellar reconnection sheets (Pen & Levin 2014). This 
complicated structure produces time and frequency vari¬ 
ations on a wide range of scales. Of course, we are con¬ 
sidering very long baselines, where the scintillation-arc 
patterns should be completely uncorrelated between an¬ 
tennas. This may lead to blurring, resulting in a 2-scale 
correlation function without particularly strong struc¬ 
ture corre sponding to the dis crete arcs seen on shorter 
baselines (Brisken et al. 2010). We do not see any direct 
evidence of scintillation arcs, as such. The magnitude of 
the visibility shows a featureless decline with increase of 
either of the 2 dimensions |r| and |/|. The GB autocor¬ 
relation functions do not show scintillation arcs either, 
for our observations. 


7. SUMMARY 

We made VLBI observations of PSR B0329+54 with 
RadioAstron at 324 MHz on projected baselines of up to 
235,000 km. Our goal was to investigate scattering by 
the interstellar medium. These properties affect radio 
observations of all celestial sources. While the results of 
such observations are in general influenced by the convo¬ 
lution of source structure with the scattering processes, 
pulsars are virtually point-like sources and signatures in 
the observational results can be directly related to the 
scattering properties of the interstellar medium. 

On long baselines, in the domain of delay r and fringe- 
rate /, the correlation function of visibility V(r,f) is a 
collection of narrow spikes, located within a region de¬ 
fined by the inverses of the scintillation bandwidth At ss 
1 / 27 tAu sc and the scintillation timescale A / ss l/ 27 rf sc . 
For shorter baselines, a sharp spike at the center of this 
region represents the average visibility; on long baselines 
where the average visibility drops to near zero, this spike 
is absent. 

The mean square visibility, ( C(t))s = (|K(r)| 2 )s, is 
well fit with a smooth model, indicating that the vis¬ 
ibility spikes are the result of random interference of 
many scattered rays. To form a quantity less subject 
to effects of noise, we convolve the mean-square left- and 
right-circular polarized visibility to form (Arl(At))s = 
(|Tl(t)| 2 |Ur(t)| 2 )s. The average correlation 

function (Art (At)) s shows two exponentials with dif¬ 
ferent characteristic timescales. The forms of ( C(r)) and 
(A"rr(t)) are well fit with a simple model, that assumes 
that the average pulse-broadening function G is the sum 
of two exponentials with different timescales. 

On 2012 Nov 26 to 29, the shorter timescale was 
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4.1 ± 0.3 /xs, and the longer timescale was 23 ± 3/xs, with 
the longer-scale exponential containing approximately 
0.38 times the power of the shorter-scale exponential. 
This double exponential may arise from anisotropic scat¬ 
tering; or from scattered radiation at large angle, perhaps 
corresponding to the subimages seen in single-dish and 
shorter-baseline observations. Further investigation of 
the properties of the image of the scattered pulsar on 
long and short baselines, using these data, will help to 
clarify the origin of the two scales. 
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ments which helped to improve the manuscript. 
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APPENDIX 

A. IMPULSE-RESPONSE FUNCTION AND VISIBILITY 
A.l. Introduction 

Under general assumptions, refraction and scattering convolve the electric held of a source with an impulse-response 
function g. This function varies with position in the observer plane, decorrelating over some lateral scale; and with 
time, as the line of sight to the source moves with respect to the scattering material, and as the scattering material 
evolves. The task of this section is to relate the impulse-response function to the statistics of visibility, as given by the 
functions C and K introduced in Section [2] above. 


A. 1.1. Notation 

The visibility Vab is the conjugated product of electric fields at two antennas (Equation [l]). We usually omit the 
subscripts indicating baseline AB on V, unless they are important for the immediate argument. We denote the Fourier 
transform from the time or delay domain (t or r) to the frequency or fringe rate domain (v or /) by and its inverse 
by 5 -1 . We accent symbols with tilde to denote quantities that depend on observing frequency v, and the same 
symbols without accent for the Fourier-conjugate domain of delay r or time t e . We assume that the variables describing 
time and frequency t e , v, r, t, /, are discrete. They range from —N/2 to N/2 — 1, where N is the number of samples in 
the time or frequency span. For r and t e one sample is the inverse of the Nyquist rate, and they can span the time to 
accumulate a single realization of the spectrum; for t and its Fourier conjugate / one sample is the averaging time for 
one spectrum, and they can span one observation. 

Our convention for normalization of the Fourier transform is that a function h{r) normalized to unit area in the 
delay domain has value unity at zero frequency: h(v = 0) = 1. Conversely, if h(v) is normalized to unit area in the 
frequency domain, h(r = 0) = 1/N. This is the “{1,-1}” convention of Wolfram Mathematica (Weisstein 2014). 
With this convention, Parseval’s Theorem takes the form: 


N/2-1 


N/2 — 1 


Y h{r)h*{T) = — Y K v )h*{v) 


t=-N/2 


(14) 


V——N/2 


A.2. Impulse-Response Function 

As noted above, the observed electric field of a pulsar E 0 b s (t 0 ) is the convolution of the electric field emitted at the 
source with a kernel g that depends on scattering: 


E 0 bs (to) 


few X T sc 

Y d(te)E s . 

t e —0 


- (to t e ) Q ® E s . 


(15) 


where we introduce the symbol (g> for convolution. The kernel g is the impulse-response function; in other words, if 
the pulsar emits a sharp spike, then the observed electric field of the pulse is simply a copy of g. Because of this 
convolution, q is al so known as the propagation kernel; it is also known as the Green’s function, and the S matrix 
(Gwinn & Johnson|[2011 and references therein). Both E src and g vary at the Nyquist rate: the inverse of the total 
observed bandwidth. Usually, we assume that the intrinsic electric field of the so urce is white n oise at the Nyquist rate: 
it is drawn from a Gaussian distribution in the complex plane at each instant (Rickett 19751. The impulse-response 
function extends over a time span of a few times t sc , representing the time over which a sharp pulse at the source 
would be received. It is zero outside this relatively narrow time window. 

If the statistics of the scattering material are stationary, the characteristic shape and scales of g will remain fixed, 
while details of amplitude and phase vary on the timescale t sc (Equation [3|. An average of the squared electric field 
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over many impulses emitted by the source over times longer than t sc will reveal the characteristic form. One simple 
model form for g that includes deterministic and random parts is the product of a non-varying envelope gD(t e ), and 
a random function g r, that varies rapidly with t e during the course of each pulse: 

g(t e ) = gD<,t e ) ■ g R {t e ,t) (16) 


Both gn and go span a few times r sc . Over that time span, g r varies wildly, and randomly; however, it exhibits nearly 
the same form for the next pulse. Over the longer timescale t sc r sc , the form of the random function gR changes 
slowly. For typical observations of a pulsar, such as those described in this paper, r sc is some fraction of the width 
of one pulse, or a few microseconds; whereas t sc is many pulsar periods, or many seconds. Such situations , where a 
convolution may have a slowly-varying kernel, are commonly treated as “dynamic spectra” (see Bracewell 2000 Ch. 

19 ). 

The intensity received by an observer for an electric-field impulse at the source, averaged over many such impulses 
with different realizations of the scattering material, is the square modulus of the deterministic part of g , which we 
call G: 


(l oba (te)) s = G(te ) = g D (te) ' 9 d(U) ( 17 ) 

Here, the subscripted angular brackets {...)s indicate a statistical average over realizations of the scattering medium, 
for example as approximated by an average over pulses spanning a time greater than t sc . Often, G is called the 
pulse-broadening function. So that propagation kernel leaves the intensity of the source unchanged, when averaged 
over time, we set: 


"22 G(t e ) = 1 


(18) 


For strong scattering, as is observed for most pulsars at most wavelengths, we expect that many different paths, with 
random amplitude and phase, will contribute to the received pulse at each instant t e . Therefore, we expect that the 
random part gR will have the statistics of a random walk at each instant: the observed electric field will be drawn from 
a circularly Gaussian distribution in the complex plane, with zero mean. On the other hand, the deterministic part go 
sets the standard deviation of g; it reflects how many paths, and with what strength, contribute at each delay. This 
model for scintillation is closely related to the amplitude-modulated-noise (AMN) model for pulsar emission (Rickett 
1975). In this model, the electric field emitted by the pulsar is the product of noise, drawn from a zero-mean Gaussian 
distribution in the complex plane, with a more-slowly varying envelope that determines the standard deviation of the 
noise at each instant. 

We suppose that the random part of the propagation kernel is completely uncorrelated in time, at the Nyquist rate, 
within its span of a few r sc . Then, at a location “A,” 


{gB,A{ie,t)g*RA^e + T A)) s — 


1 if r = 0 
0 if t 0 


(19) 


On the other hand, gRA is nearly the same for each emitted pulse; it changes only over the longer timescale t sc . The 
question of how this slower variation of gn with time depends upon baseline length is much more complicated, and we 
discuss it briefly below. However, if the lateral separation of the two stations A and B is much greater than the scale 
of the scattering pattern, then the random parts of g for the two stations, gR A and 9rBi are completely uncorrelated. 


A.3. Visibility: Dynamic Cross-Power Spectrum 

As the previous discussion shows, the impulse-response function involves three timescales: gn changes at the Nyquist 
rate; go varies over the typical span of the impulse-response function r sc ; and the time for the random variatio ns of qn to 
chang e is t sc . The dynamic cross-power spectrum provides a useful description for these different variations (Bracewell 
2000). A single sample of the cross-power spectrum, when averaged over time less than t sc , has the characteristic scale 
Ar' sc ss 1/2itt sc , resulting from the finite span of the impulse-response function and the uncertainty principle. The 
time variation of the cross-power spectrum over times t > t sc captures the changes of g r. 

Visibility V in the domain of frequency v and time t is the product of the Fourier transforms of electric fields at 
stations A and B (see Equation [I]): 


V A B{v,t) = E A {v,t)E* B {v,t) 


( 20 ) 


We suppose that each sample of the cross-spectrum is averaged over many realizations of the source electric field E src , 
over a time short compared with t sc . This reduces noise from the source and backgrounds. 

One may represent the visibility in four domains, linked by Fourier transforms of frequency v to delay r, and time 
t to fringe rate /, as Figure [I] illustrates. In this paper we are particularly concerned with Vab(t, t), visibility in 
the domain of delay r and time t. This is the Fourier transform of Vab(v, t). The convolution theorem for Fourier 
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transforms shows that V (r, t ) is the cross-correlation function of electric fields in the time domain: 


N/ 2-1 

Vab{ t A) = 'y ] E A (t + te)Eg(t + t e + r) 

t e —-N/2 


(21) 


Here, t e indexes individual samples of electric field, over a short interval near the index time of the measurement of 
the cross-power spectrum, t. 

Visibility in the delay-rate domain, conjugate to the frequency-time domain, takes the form: 


V(tJ) = 5 1 


V — 


3t- 


V(v,t) 


( 22 ) 


Searches for interference fringes are often conducted in this domain: because absolute calibration of delay and fringe 
rate are usually impossible for very-long baseline interferometry, the peak of I V(t, f) I can be used to determine them 
(Thompson et al. ]§007). 

A.3.1. Visibility and Impulse-Response Function 

Visibilit y de pen ds o n the separation b of stations A and 1?, as well as on delay and rate, or time and frequency. From 
Equations |15| and |21[ and the assumption that the electric field of the source E src is a stationary random variable 
without correlation in time, we find that visibility in the delay domain is the cross-correlation function of g at the two 
stations: 


V A b(t) = g A (n) <8> ffs(-«) 

K,—YT 


(23) 


This leads to the expected form of V{r): a spike at r = 0, with average magnitude equal to the average correlation 
Pab\ and a broad component of width t sc , with random amplitude, and phase variations that increase with baseline 
length, corresponding to the random character of gn and its decorrelation with increasing baseline. Equations [l9| and 
[Tsl show that pab = 1 for pra = Prb, and pab —> 0 for uncorrelated qra and gas- Figures [2] and [3] show examples. 

On intermediate baselines, the time structure of the correlation of gn is more complicated) in a way that depends 
on the geometrical distribution of the paths that contribute to gR. For scattering material concentrated in a thin 
screen, for example, the shortest-length paths result in small t e in the impulse-response function, and also tend to 
appear at small angles at the observer. Thus, at small delays correlation is high even for rather long baselines; whereas 
at long delays correlation is poor even for shorter baselines. Thus, correlation between antennas should decrease at 
later times t e within gR. This correlation is imprecise for scattering material distributed along the line of sight, where 
many deflections along the line of sight lead to a large time lag t e , but little or no angular deflection at the observer’s 
interferometer. Moreover, in the frequency domain, dynamic single-dish spectra can show slants and complicated 


patterns (Hewish 


1980 Stinebring et al. 2001), suggesting complicated correlations of time and delay in the observer 


plane. Anal ysis~ of t he v isibility is thus easiest on very short baselines and very long ones. 

Equations [22] and [23] provide the relation of the visibility in the frequency-time domain to the impulse-response 
function: 

Vab = (go ® 9ra) ( g* D ® g*R B ) (24) 

Averaging the visibility over many scintillations yields the average visibility, pab- 

Vab^J s = Pab (25) 

A.3.2. Square Modulus of Visibility C and Correlation Function K 

The average of V (r) over many realizations of scintillations leaves the delta-function at the origin that corresponds 
to the average visibility: 


(Vab(t)) s = 0 + pab6° t (26) 

Here, Sq is the Kronecker delta-function, with value 1 if r = 0 and 0 otherwise. 

The secondary spectrum C(t, /) may be defined as the square modulus of V (r, /a 

C(t,/) = = |F(r,/)| 2 (27) 

This function provides information similar to \V(r, f)\, shown in Figure [2j but is easier to deal with statistically. 
Because C(r, /) is the square modulus of the complex visibility V(t , /), background noise adds noise to C, except for 


9 |Brisken et al.j |2010 ) define the secondary spectrum as 
V (t, f)V(—r, —/). For zero baselines V(v,t) is real, so that 
V(t, f) = V*(—r, — /), and our and their expressions are identical. 
Their expression includes phase information in an elegant way for 
their short baseline, where departures from zero phase are small. 
For observations on long baselines, their expression is impractical 


because identification of the origin of (r, /) is not possible, as Fig- 
ure[2]shows. Consequently, the pair (r, /) and (—r, — /) cannot be 
combined reliably. They also use the accent ~ for visibility in the 
delay domain, although not for C; whereas we use the accent for 
quantities in the frequency domain of u. 
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the central lag, where it may contribute a constant offset. Similarly, self-noise will add noise with an envelope that 
follows the average form of C. This behavior is in contrast to that of correlation functions of single-dish quantities 
such as the intensity, where noise can contribute to the mean correlation function. 

Our long baselines fully resolve the scattered image. The phases of scintillation elements in V(is, t) appear to be 
random, and the phases of C(t, /) show no discernible patterns. An inverse Fourier transform from fringe rate / to 
time t leads to C(r,t): 

C(T,t) = r f \ t [C(T,f)}. (28) 

Evaluated at the fringe rate fmax of its peak magnitude, the secondary spectrum C(t, f ma x) is a time average, that 
approximates an average over realizations of the scintillation pattern of C(r,t ): 

(c(T)) s aC(r,U) (29) 

The autocorrelation function of C is K: 

K(At) = C(t, t)C(r + At, t ) (30) 

Note that C(r,t ) appears without averaging in this expression. Conveniently, the correlation Krl(At, t) between the 
secondary spect ra C i n the right and left circular polarizations eliminates some effects of noise and interference, as 
noted in Section [5.2.31 

The behaviors of C and K are simplest to describe on very short baselines, where 3 a(n) = ffs(ft) for all k, and 
Pab = 1; and for very long baselines, where the random parts of ^(k) and 3 b(k) are completely uncorrelated, so that 
Pab = 0. If the correlation of the random parts of the propagation kernels pab is constant, then the ensemble-average 
values of these correlation functions are: 

(c(T)) s = c *{av G ) + (p 2 Alj )d; ] (31) 

(k( At)) s = ^ (G ® G_ ® G ® G_) (1 + 4 p\ B ) + 2 + p\ B ^j 5^ (32) 

Again, is the Kronecker delta-function, with value 1 if r = 0 and 0 otherwise; and similarly for At. The time- 
reversed pulse-broadening function is G_, given by G-(k ) = G(—k). The constant Co is the mean square of G: 

= E t(: G 2 . 


B. TWO EXPONENTIALS 


We consider a situation where the impulse-response function is the sum of two exponentials with different time 
constants k\ and k 2 , with a rapid rise from G = 0 at t = 0. As we discuss below, this may result in a variety of 
circumstances. We parametrize the impulse-response function: 


f Ai/cie- felT + A 2 k 2 e~ k2T r> 0 
\0 t < 0 


(33) 


In the text, we also make use of the inverse scales 7i = 1/Aq and t 2 = l//c2; these can provide better physical insight. 
The autocorrelation of G provides the form for ( C(r))g , the mean square visibility in the delay-time domain, as given 
by Equation |31 1 For this impulse-response function, under the assumption that the baseline is so long that pab = 0, 
this takes theform: 


^(c(Ar)) 

C 2 \ / s \^i + k 2 


M 

2 


a -fci|Ar| 


A 2 k 2 


Aiki 

k\ + k 2 


A 2 

T 


0 -fc 2 |Ar| 


(34) 


The autocorrelation function K of the mean square visibility for this impulse-response function then takes the form: 

N 

C 2 \ " -is 

+ (((or — P)ki + ( a 2 + (3)k 2 ) 5^ 


(a(|At|)) = (aifci(l + fci |Ar|) - pk 2 ) e _fe il A ' r l + ( a2 k 2 (i + k 2 |Ar|) + ph) e ~ k ^ AT ' 

: ' ' S 


(35) 


where again S At is the Kronecker delta function, and: 


Oil = A\ 
a 2 = Al 



A 2 k 2 \ 
k\ + k 2 J 

Aiki y 
k\ + k 2 ) 


AiA 2 kik 2 (2A 2 k 2 + A\(ki + k 2 ^(2Aiki + A2(Aq + k 2 )) 


(36) 


(h - fc 2 )(fci + k 2 ) 3 
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Thus, the two exponential scales l/ki, l/k 2 appear again, in K , as do the weights A±, A 2 . 

C. ANISOTROPIC SCATTERING IN A THIN SCREEN 

An observer sees an anisotropic distribution of radiation from a screen at distance D from the observer. The 
probability of receiving radiation from the screen at position ( 0 X , 0 y ) is: 


P(0x,0y) = 


1 


‘Z'K CT x CTy 


exp ■ 


dd x ddy 


(37) 


We suppose without loss of generality that a y > <j x . If the source is at infinite distance beyond the screen, the delay 
t along this path (ignoring any contribution from the screen) is: 


ct = d (e 2 x + e 2 y ) 


(38) 


where c is the speed of light. If the source is at distance R beyond the screen, then D is replaced by RD/(R + D) 
in this and subsequent equations. We convert the distribution of angles in Equation 37 to coordinates {r,<p), where 
(j) = arctan (0 y /6 x ). The resulting distribution of (r, <j>) is: 


P(t, 4>) d(j> dr = y/l + a 2 ^^2 ex P j - 0]Ja 2 + cos2 | ^ 

where a 2 = a 2 /a^ — 1 parametrizes the anisotropy. We integrate over <f> to find the pulse-broadening function: 


G{t) dr = P(t, cj)) 
Jo 


d(j)dr 


Thus, 


G(t)cIt = \/l + a 2 exp / 0 ( 


In 


2 r 2 J t 2 


(39) 


(40) 


(41) 


Here, Iq is the regular modified Bessel function of order 0, and CT 2 = 2 Da 2 . Note that for a —> 0, this distribution 
becomes the familiar exponential distribution, with scale r 2 , as expected for a x = a y . Rickett et al. (2009) present a 
similar expression. 

In the case of a > 0, at small r the distribution has the limit: 


lim G(t) dr^V l + a 2 exp (- (2 + Q2) -l — 
r-y 0 [ 2 r 2 J r 2 


At large values, 


so that 


lim Iq(u) 

U—> OO 


27 TU 


lim G(t) dr —► \J\ + a 2 * 


1 


T2 


2tt(2 + a 2 ) t 


dr 

exp i-> — 

T 2 J t 2 


(42) 


(43) 


(44) 


At a particular scale r, the logarithmic derivative of G(t ) is I/ 72 , although the coefficient depends weakly on r. Thus, 
G(t) exhibits two exponential scales: t 2 at large r, and t\ = t 2 /{\ + a 2 /2) at small r. The relative strength of the 
scales is about 


Ai 

A 2 


1 + a 2 


27t(2 + a 2 ) 

near the larger scale t 2 , with the larger scale being the weaker. 


(45) 


REFERENCES 


Andrianov, A. S., Guirin, I. A., Jarov, V. E., et al. 2014, Vestnik 
NPO, 24, 55 

Bartel, N., Ratner, M. I., Shapiro, I. I., et al. 1985, AJ, 90, 2532 
Boldyrev, S., & Gwinn, C. R. 2003, Physical Review Letters, 91, 
131101 

Bracewell, R. N. 2000, The Fourier transform and its applications 
(McGraw Hill) 


Brisken, W. F., Benson, J. M., Goss, W. M., & Thorsett, S. E. 
2002, ApJ, 571, 906 

Brisken, W. F., Macquart, J.-P., Gao, J. J., et al. 2010, ApJ, 708, 
232 

Britton, M. C., Gwinn, C. R., & Ojeda, M. J. 1998, ApJ, 501, 
L101 

Desai, K. M., Gwinn, C. R., Reynolds, J., et al. 1992, ApJ, 393, 
L75 




















RadioAstron discovery of scaterring substructure in PSR 0329+54 


15 


Edwards, R. T., Hobbs, G. B., Sz Manchester, R. N. 2006, 
MNRAS, 372, 1549 

Goodman. J., Sz Narayan, R. 1989, MNRAS, 238, 995 

Gwinn, C. R., 2001, ApJ, 554,1197 

Gwinn, C. R., Hirano, C., Sz Boldyrev, S. 2006, ASzA, 453, 595 

Gwinn, C. R., Sz Johnson, M. D. 2011, ApJ, 733, 51 

Gwinn, C. R., Johnson, M. D., Smirnova, T. V., Sz Stinebring, 

D. R. 2011, ApJ, 733, 52 

Gwinn, C. R., Kovalev, Y. Y., Johnson, M. D., Sz Soglasnov, 

V. A. 2014, ApJ, 794, L14 

Gwinn, C. R., Johnson, M. D., Reynolds, J. E., et al. 2012, ApJ, 
758, 6 

Hagfors, T. 1976, Methods of Experimental Physics, 12, 119 

Hewish, A. 1980, MNRAS, 192, 799 

Johnson, M. D., Kovalev, Y. Y., Gwinn, C. R., et al. 2016, ApJ, 
820, L10 

Kardashev, N. S., Khartov, V. V., Abramov, V. V., et al. 2013, 
Astronomy Reports, 57, 153 

Kondratiev, V. I., Popov, M. V., Soglasnov, V. A., et al. 2007, 
Astronomical and Astrophysical Transactions, 26, 585 

Lorimer, D. R., Yates, J. A., Lyne, A. G., Sz Gould, D. M. 1995, 
MNRAS, 273, 411 

Narayan, R. Sz Goodman, J. 1989, Philosophical Transactions of 
the Royal Society of London Series A, 341, 151 

Narayan, R., Sz Goodman, J. 1989, MNRAS, 238, 963 


Pen, U.-L., & Levin, Y. 2014, MNRAS, 442, 3338 

Popov, M. V., Sz Soglasnov, V. A. 1984, Soviet Ast., 28, 424 

Prokhorov, A., Bunkin, V. F., Gochelashvily, K. S., Sz Shishov, 
V. I. 1975, Proc. IEEE, 63, 790 

Rickett, B. J. 1975, ApJ, 197, 185 

—. 1977, ARA&A, 15, 479 

Rickett, B., Johnston, S., Tomlinson, T., Sz Reynolds, J. 2009, 
MNRAS, 395, 1391 

Semenkov, K. V., Soglasnov, V. A., & Popov, M. V. 2004, 
Astronomy Reports, 48, 457 

Shishov, V. I., Smirnova, T. V., Sieber, W., et al. 2003, A&A, 
404, 557 

Smirnova, T. V., Shishov, V. I., Popov, M. V., et al. 2014, ApJ, 
786, 115 

Stinebring, D. R., McLaughlin, M. A., Cordes, J. M., et al. 2001, 
ApJ, 549, L97 

Thompson, A. R., Moran, J. M., Sz Swenson, G. W. 2007, 

Interferometry and Synthesis in Radio Astronomy, John Wiley 
& Sons, 2007. (Wiley) 

Weisstein, E. W. 2014, Fourier Transform, 

http://mathworld.wolfram.com/FourierTransform.html 
accessed: ^Ui0-Uy-Uli 

Yangalov, A. K., Popov, M. V., Soglasnov, V. A., et al. 2001, 
Ap&SS, 278, 39 


